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Abstract 

We calculate the differential cross section for semi-inclusive pion production in 
electron proton reactions using a model where the physical quark fluctuates with 
some probability to quark plus pion. The kinematic regions for a determination 
of this ‘pion cloud’ are evaluated. 


1. Introduction 

Since Yukawa the nucleon has been perceived as an extended object with a pionic 
cloud around it. The one-pion exchange interaction between nucleons is a natural 
consequence of pion emission and absorption by nucleons. The extent of the pionic 
cloud of the nucleon has, however, been debated. Proponents of the N-N-one-boson 
exchange model use form factors with cut-off values A ss 2 — 2.5 GeV, which give a 
large amount of pion admixture to the nucleon. In an analysis of pion-electroproduction 
data on the proton [1] an estimate of the pion cloud of the nucleon has yielded a rather 
small value, namely a 3% admixture of n7r + to the naked p state. This small probability 
would make it unrealistic to use coherent pion emission of the nucleon to create sea 
quarks in the nucleon, which is a process proposed in Refs. [2], 

On the other hand, it is very natural to couple quarks to scalar (<r) and pseudoscalar 
mesons (if) in a chirally invariant way. The first coupling leads to the constituent quark 
mass after spontaneous breaking of the (<r, if )-symmetry. The second coupling to pions 
produces a pion cloud of the constituent quark. There is a difference between a pion 
cloud of the n ucleon and the pion cloud of the quark. Emission and absorption of a pion 
by the same quark determines the pion cloud of the constituent quark. The interme¬ 
diate states on the nucleon level would be | 'iqjr) states of arbitrary excitation energy. 

1 Present address: Alfred Weber- Institut, Grabengasse 14, D-69117 Heidelberg 
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Allowed intermediate states of the nucleonic pion cloud are | N*tt) or | Att) states. In 
this sense the pion is part of the hadronic system. The analysis of electroproduction at 
small inclusive mass and correspondingly z = E^/v —> 1 will give only the pion cloud 
of the nucleon, as described in ref. [1], We suggest that the constituent quark structure 
is the more interesting quantity to study, as was also proposed in Ref. [3]. 

The constituent quark picture has been a simple and convincing concept to explain 
the static properties of baryons and mesons. It is surprising how well magnetic moments 
and masses of baryons are reproduced. In deep inelastic lepton nucleon scattering the 
quarks are visible with their momentum distributions and charges. Yet the constituent 
quark model has never been able to explain the measured structure functions of the 
proton. Especially the transition with Q 2 from an extended constituent quark with its 
own internal structure to a pointlike current quark has not been investigated. 

At low Q 2 , the constituent quark has intrinsic structure, namely a pion and possibly 
also a gluon cloud. The baryon will be described by a light cone wave function of three 
constituent quarks. This wave function takes care of the relativistic motion of the 
quarks and can, in principle, be obtained from a suitable light cone Hamiltonian. 

In this article, we want to focus on the phenomenological aspects of the constituent 
quark. We want to define a procedure to separate its internal structure from the bound 
state dynamics contained in the quark wave function. 

For this purpose, we first analyze in Section 2 the quasi-static properties of the 
nucleon which are derived from moments of deep inelastic structure functions, namely 
the Gottfried sum rule, Sg, the axial vector coupling, and the spin sum rules, T^’A 

These moments depend on integrated pion probability distributions. In the second 
step (Section 3) we determine the semi-inclusive pion production cross section which 
gives differential information about the pion cloud of the constituent quark. Section 
4 analyzes the experimental possibilities to measure the semi-inclusive process in the 
relevant Q 2 , x B , and z range described. Finally we conclude in Section 5. 

2. The Pion Content of the Constituent Quark 

In the naive parton model the proton consists of two up-quarks and one down-quark: 
|p) = | uud). Beside the valence quarks there is a sea of quark antiquark excitations 
which, however, does not affect the sum rules concerning electric charge, baryon number 
or strangeness. In an obvious notation, the total content of flavor i in the proton can 
be written as 

qi(x) = q t , v (x) + qi, s (x) + qi, s (x) where i = u,d,s . (1) 
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Making two assumptions, namely 

1. Isospin invariance: 

u p (x ) = d n (x ) and d p (x ) = u n (x ) (2) 

2. Flavor symmetry of the sea: 

u s (x) = u s (x) = d s (x) = d s (x) (3) 

the Gottfried sum rule [4] is derived in this simple model to be 



o 

Measurements at CERN [5] have yielded an experimental value of 

Sq P = 0.240 ± 0.016 (5) 

which implies a substantial violation of the Gottfried sum rule, Eq. 4. This so-called 
GSR defect can be accommodated in two ways. It has been suggested [7] that contri¬ 
butions from the unmeasured region of small x (x < 0.004) could make up this defect. 
However, recent measurements of the structure functions for very small x (down to 
x ~ 10 -4 ) at HERA [8] show that the NMC extrapolation employed in Eq. 5 is fairly 
good. The other interpretation, first proposed by Eichten, Hinchliffe and Quigg [9] and 
assumed in this paper, is that assumption 2 used to derive the Gottfried sum rule is 
wrong and the sea of light quarks is flavor asymmetric instead, i.e. 

u s (x) = u s (x) ^ d s (x) = d s (x) . (6) 

Now we get a prediction for the Gottfried sum which is different from the value 1/3 
depending on flavor asymmetry: 

1 2 r 

Ag = - + - y dx \u s (x) - d s (x)^ . (7) 

o 

The observed value of Sg = 0.240 ± 0.016 means that in the sea of the proton there 
is an excess of d-quarks over u-quarks [10]. The origin of this flavor asymmetry can 
be understood in the framework of effective chiral quark theory [11,12], Quarks (that 
is: massive constituent quarks) can fluctuate into quarks and pions, the only bosonic 
degrees of freedom considered here. Hence, the physical constituent quark is the co¬ 
herent sum of a single naked (albeit massive) constituent quark and a state in which 
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the quark is dissociated into constituent quark and pion. It is this idea that we want 
to term “pion cloud of the constituent quark”. 

Before presenting a more detailed discussion in the next section of the properties of 
this pion cloud of the constituent quark, we want to employ a toy model [9] that already 
captures the essential features. In this simplified picture a nucleon consists of up- and 
down-quarks only. In the chiral quark model a u-quark can emit a tt + (containing a 
valence u-quark and a valence d-antiquark) or a 7r° as depicted in Fig. 1. The model 
pa ram el ei! describing this idea is 

a = | ( d7r + 1 u ) | 2 , (8) 

the probability for an up-quark to turn into a down-quark with the emission of a 7r + . 
We assume that this fluctuation is small enough to be treated as a perturbation. 



Figure 1. Fluctuation of a u-quark into quark and pion. 


In this picture a physical constituent quark (U) after one such interaction is a 
superposition of different (naked) states: 
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The flavor content of the proton and neutron in this simplified model hence is given by 
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( 12 ) 

One sees: by assuming a pion cloud of the u- and d-constituent quarks (that is: a ^ 0) 
there arises a flavor asymmetry in the nucleon’s sea the extent of which depends on 
the value of the model parameter a. Hence, the GSR defect can be expressed in terms 
of a: 


1 2 a 




3 3 


(13) 


Of course, the flavor asymmetry represented by the parameter a also affects other static 
and quasi-static properties of the nucleon. Let 

l 

Aq t := J dx \q](x) - ^(x)] (14) 

o 


denote the integrated difference of quarks (and antiquarks; see Eq. 1) with helicity 
parallel (|) and antiparallel (J,) to the nucleon’s helicity. In the naive quark-parton 
model with SU(6)-wave functions for a polarized proton one has 


A u 
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4 
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( 15 ) 


and gets the following prediction for the axial vector coupling constant gA 


gA = A u — Ad 


5 

3 ' 


( 16 ) 


If we consider a pion cloud of the constituent quark in the form of the just explained toy 
model, the spin content of the constituent quark is changed, since the pion is emitted 
in a p -wave state. One gets 
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(17) 
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and thus a value of g A 


g A = AU-AD = ^(l-a) , (18) 

which is modified as compared to the naive expectations in Eq. 16 by the existence of 
a pionic admixture (a ^ 0) to the constituent quark. 

Similarly, g A , defined as the proton matrix element of the flavor singlet axial current 

(p\ V’7 75'0 \p) =■ 9aX ] 0X, ( 19 ) 


where ^ = (u,d,s) denotes the quark held and | p ) is the wave function of a proton at 
rest with Pauli spinor y, can be expressed using Eq. 17: 

g° A = AU + AD + AS = 1 - 3a . (20) 


g A measures how much of the proton’s helicity is carried by constituent quarks. In 
the naive quark parton model a = 0 and the proton’s helicity is carried entirely by 
quarks (g A = 1). In the model of a pion cloud this value is reduced by the pions’ (or 
sea-quarks’) contribution and thus g A < 1. 

A rich source of insight into nucleon structure are structure functions for deep 
inelastic inclusive scattering on polarized nucleon targets. They can also be expressed 
using quark distribution functions: 

9i(x) = [^( x ) ~ ^( x )\ • ( 21 ) 

i 


For experimental reasons, instead of the structure functions one rather considers mo¬ 
ments of them: 

l 

Ti := J dxg^x) . (22) 

o 

In the pion cloud model one obtains for the proton and neutron spin structure functions 
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(23) 

(24) 


We will now analyse the available experimental data on the above mentioned quantities 
Sg, 9a-, 9 a and to determine the value of a. 

The Gottfried sum has been measured by the New Muon Collaboration (NMC) at 
CERN [5]. They determined the ratio F^ 11 /in the kinematic range 0.004 < x < 0.8 
and 0.4 GeV 2 < Q 2 < 190 GeV 2 and obtained 


5g(0.004, 0.8; (Q 2 ) = 4 GeV 2 ) = 0.227 ± 0.006 (25) 
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This value has been extrapolated to x —> 1 and x —> 0 assuming Regge behavior for 
small x with the result 

S G = 0.240 ± 0.016 . (26) 

An improved analysis of the data [6] has yielded the value 

S G = 0.235 ± 0.026 (27) 

which we will employ in our analysis. This improved data analysis also suggested that 
a so-called shadowing effect occurs when scattering on deuterium targets: the neutron 
structure function at small x is larger than the difference between the deuterium and 
the proton structure functions. This effect has not yet been measured directly, but has 
to be estimated theoretically. Since NMC obtained the neutron structure function as 
the difference between deuterium and proton structure functions its result for S G has 
to be corrected for this shadowing effect. Various calculations on the shadowing effect 
have already been done [13-15]. For the measured range x > 0.004 the authors of [13] 
obtain a correction to S G of fr shadow g G = —0.026. For the whole range 0 < x < 1 the 
authors of [14] give a correction between —0.010 and —0.026. However, comparing these 
calculations with data for photoproduction shows that the effect is underestimated in 
[13]. Using a different model for diffractive scattering [15] one obtains a correction to 
S G in the measured range x > 0.004 of fr shadow g G = _0.043 and an additional correction 
in the region not measured (x < 0.004) of $ shadow S G = _0.038. For further analysis we 
employ the latter result and thus have 

A g ( 0, 1; (Q 2 ) = 4 GeV 2 ; corrected) = 0.154 ± 0.026 . (28) 

Note that the correction for shadowing amounts to 30 % of the original NMC result! 

The quantity known with by far the highest precision [16] is the axial vector coupling 
constant 

g A = 1-2573 ± 0.0028 . (29) 

The spin structure functions gi(x) of the proton and the neutron have been mea¬ 
sured recently by several collaborations. The European Muon Collaboration (EMC) 
has measured the moment of the proton structure function in the range Q 2 > 10 GeV 2 
(with ( Q 2 ) = 10.7 GeV 2 ) by scattering muons off hydrogen targets [17]: 

Tl = 0.126 ± 0.025 . (30) 

At CERN the Spin Muon Collaboration (SMC) [18] has measured the moment of the 
deuteron structure function in the kinematic range 0.006 < x < 0.6 and 1 GeV 2 < 
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Q 2 < 30 GeV 2 (with ( Q 2 ) = 4.6 GeV 2 ) to be 


= 0.023 ± 0.035 . (31) 

At SLAC the E142 experiment has measured [19] the spin structure function of the 
neutron in the kinematic range 0.03 < x < 0.6 at ( Q 2 ) = 2 GeV 2 to be 

0.6 

f dxg^(x) = -0.019 ±0.012 . (32) 

0.03 

From this value an extrapolation yields 

l 

17 = J dx g?(x) = -0.022 ± 0.011 . (33) 

0 

Ellis and Karliner [20] pointed out that extreme care should be taken when comparing 
these results on polarized scattering because the experiments have been carried out at 
different Q 2 . They performed a careful analysis of the experimental results and ap¬ 
plied various corrections ((^-dependence, perturbative QCD corrections, higher twist 
effects). They find that the data from E142 and SMC are consistent within less than 
one standard deviation with the Bjorken sum rule (/ d x (gl(x) — g™(x)). One sees 
that the (^-corrections are essential for understanding the data. For this reason we 
will employ the corrected values to determine a. Using all available data on polar¬ 
ized lepton-nucleon scattering and applying the above mentioned corrections Ellis and 
Karliner obtain a best fit to the quark contribution to the nucleon helicity: 

AS = g° A = 0.37 ± 0.07 . (34) 

Table 1 summarizes the experimental results on the various quantities, their relation 
with the parameter a and the result for a which is obtained from the respective values. 

Fig. 2 shows the various values for the model parameter a. We also plot as a 
thick vertical line the mean from the results of the three independent experimental 
determinations, Eqs. 28, 29, and 34. The best known quantity is g^- It is this value 
which therefore essentially determines the weighted mean for a: 

a = 0.246 ± 0.002 . (35) 

As one sees in Fig. 2 the various corrections are essential. By correcting the NMC 
result on Sg for shadowing effect and the E142 result for (^-dependence the original 
values are shifted to within la of the weighted mean of a = 0.246. We take this as 
evidence that the employed toy model captures at least the essential features of the 
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Observable 

Theory 

Measurement 

Result for a 

S G (NMC) 

A g = 

(l/3)(l-2o) 

0.235 ±0.026 

0.148 ±0.039 

S G ± ± hadow ± G 



0.154 ±0.026 

0.269 ±0.039 

9a 

9 A = 

(5/3)(l -a) 

1.2573 ± 0.0028 

0.246 ±0.002 

H) (EMC) 

r; = 

(5/18)(l - 2 a) 

0.126 ±0.025 

0.273 ±0.045 

r? (smc) 

rn 

1 1 “ 

—(5/18)a 

-0.08 ±0.08 

0.29 ±0.29 

r? (E142) 



- 0.022 ± 0.011 

0.079 ±0.040 

T ] 1 (E142, corr.) 



-0.056 ±0.011 

0.202 ±0.040 

best fit for g° A 

ii 

1 — 3a 

0.370 ±0.070 

0.210 ±0.023 


Table 1. Determination of a by various observables. 


idea of a pion cloud of the constituent quark and that this idea is capable of explaining 
some properties of the nucleon. 

However, it might come as a surprise that the probability for a 7 r + in an up- 
constituent quark turns out to be almost 25 %. A perturbative calculation in the 
linear cr-model with pseudoscalar quark-pion-coupling (coupling constant g vqq = 3.76) 
reproduces this value if one uses a regularization parameter of A = 3.37 GeV. This 
seems to be very large for a hadronic cut-off. A possible explanation for both large 
values might be the non-relativistic character of the employed models. It should be ap¬ 
propriate to describe the nucleon as consisting of three relativistic quarks. Relativistic 
quark models with light-cone wave functions for the constituent quarks have already 
been used for a quantitative analysis of static and quasi-static properties of the nucleon 
[21,22], While such calculations can predict static properties of the nucleon, like <^ 4 , 
almost exactly without recurring to an effective description like a pion cloud, they nev¬ 
ertheless cannot give a satisfactory agreement with experiment as far as quasi-static 
properties are concerned. 

We conclude that even in a relativistic description of what constitutes a nucleon 
the idea of a pion cloud of the constituent quark might be helpful. However, the above 
performed estimates of the pion probability a in the constituent quark, which are based 
on integral properties of the nucleon, are only of limited use. A direct and differential 
measurement of the pionic substructure of constituent quarks will offer more insight. 
To this purpose we turn in the next section to a differential description of the nucleon 
as seen in semi-inclusive electroproduction of pions on nucleon targets. 
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Figure 2. The various results for a from different observables. The thick vertical line 
indicates a = 0.246. 

3. Semi-Inclusive Pion Production 

We consider electron scattering on the proton with at least one pion 7 r(k) in the final 
state: 

e + p —>■ e' + 7 T + X . (36) 

Semi-inclusive pion production can be looked at as a soft hadronization process. 
Data [23] at high energies agree rather well with LUND-simulations. However, in 
certain kinematic regions higher twist effects can visibly modify the leading gluon cas¬ 
cade. For example, one expects contributions [24-26] from the hard fragmentation 
7 *q —> irq. 

We will describe the pion production in the naive parton model with the follow¬ 
ing assumptions. The nucleon consists of three valence constituent quarks, which are 
weakly bound. The square of their wave functions integrated over transverse momen¬ 
tum p±_ up to the resolution Q 2 gives the probability to find a quark with flavor q and 
light cone momentum fraction x in the proton: 
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The constituents are massive Dirac particles ( Mq = g vqq (<r)) with a pointlike coupling 
to the photon and a pseudoscalar coupling ig^qq ^ 1° the pion. The coupling 

constant g^qq couples the quark equally to both fields of the (|,|) SU(2)iXSU(2)# 
multiplet (<r, fr). The coupling strength is determined from the empirical value of the 
constituent quark mass Mq = 350 MeV and the value of (<r) = = 93.3 MeV: 

g^qq = 3.76 . (38) 


Equally one could have estimated g vqq from quark additivity and the 7r-N-coupling 
constant g^NN = 13.5: 


_ g^NN ., -- 

” ~ 3 9 . 4 ' 


(39) 


In the parton model, pion production is calculated as a folding of the partonic cross 
section eq —> e'q'ir with the distribution N q (x) of the quarks: 
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da ep ~ 


e' ttX 


= £/ N q (x)2E T 


dcr eq ^ e q v 


dx 


(40) 


d 3 kdx B dy q J d 3 kdx B dy 

The kinematical variables of this semi-inclusive process are on the electron side 
x B = Q 2 /2Mjvu and y = v/E and on the pion side k p = (E V} k). It is useful to discuss 
the kinematics in the next section in more detail. 


3.1. Kinematics of Semi-Inclusive Pion-Production 


In Fig. 3 we give a definition of all momenta in the target-at-rest system, P p = (Afjy, 0 ). 
In the semi-inclusive reaction there are two planes of interest, the lepton scattering 
plane (7,7') and the photon-pion plane (<f, k). They include an angle between them. 
The energy of the incoming electron is E } the energy of the outgoing electron is E'. 
The photon energy is z/, the pion energy is E v . The pion momentum k is divided into 
a momentum parallel to the direction of the photon momentum q 


k ■ q 




l?l 


and transverse to the photon momentum q 


(41) 


k L 


k — k\\ 


(42) 


Unpolarized cross sections do not depend on the azimuthal angle of the scattered 
electron, therefore two invariants are sufficient to characterize the virtual photon. These 
are 


Q 2 = ~q 2 
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4 EE' sin 2 d e /2 


P ■ q 

= E-E 
M n 


(43) 
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Figure 3. The reaction ep —> PttX in the target-rest- system. 


or 


x B = 


y = 


-q 


Q 7 


2 P ■ q 2 Afjv v 
P ■ q v 

~P~l ~ E ' 


(44) 


The momentum of the pion has three components. We separate the azimuthal angle 
from the two residual variables k\\ and \k±\. The parallel momentum k\\ can be 
replaced by the energy variable E v or the scaling variable 

P ■ k Ei r 
P ■ q v 


z = 


(45) 


A conventional set of variables thus consists of ( Q 2 , x B , z, \k]_\ 
Another invariant variable besides Q 2 , P ■ q , and P ■ k is 


q. k = v l E v - \ 1 + 
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(46) 
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It can be simplified for £ = / Q 2 <C 1 , C = + m l/Q 2 G and CC 2 2 t° 


q ■ k 


1 

2 


z x _g 



( 47 ) 


which shows that in the Bjorken limit, Q 2 —>■ oo, q ■ k = — ^Q 2 z is no longer an 
independent variable. 

It is advantageous to introduce the Mandelstam variables of the parton subprocess 
J*(q) + q(p) —> 7 r(k) + q'(p') (with q ■ k from Eq. 46): 

s = (P + q ) 2 = M| + Q 2 (--l) 

V x B / 

h = (q — p' f = Mq + Q 2 + 2q ■ k 

t = (q — k) 2 = — Q 2 + ml — 2q ■ k . (48) 


Note that p and p' denote the quarks’ 4-momenta, whereas P and P' are nucleon 
4-momenta. The Mandelstam variables are related by s + u + t = 2 Mq + to 2 — Q 2 . 

Since the incoming photon is a virtual particle with q 2 = — Q 2 < 0, the kinematical 
boundaries for the variables x B , z, Q 2 , and k\ are interesting. The condition, that the 
quark after the production is on-mass-shell, 


p ' 1 = (,T P + g - Ilf = Ml , 


(49) 


expresses the quark momentum fraction x in the nucleon to 
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(50) 


From the condition x < 1 follows a maximum x B for a given set Q 2 , z, and kf . In Fig. 4 
we show the maximal possible value of \k±\ for fixed Q 2 = 1 GeV 2 . One sees that the 
allowed region for x B becomes maximal for intermediate z values. The maximal | A;j_ | is 
obtained for small x B . 


3.2. Parametrization of Semi-Inclusive Pion Production 


The inelastic electron scattering cross section can be factorized into a leptonic and 
hadronic tensor multiplied by the photon propagator squared 
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Figure 4. Maximal value of |A?j_| as a function of ~ and x B at a fixed value of Q 2 = 
1 GeV 2 . 


where 
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and 
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W- w is parametrized by four independent structure functions T4y(,r B , Q 2 , kj _) assum¬ 
ing parity and time reversal invariance [27] 
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with the polarization vectors 
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(54) 


(55) 
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The cross section for definite photon polarizations can be obtained by contracting W^ v 
with e tl e v of the wanted sort, e.g. Wl = e n{L) W^ v e v (L) for the longitudinal photon 
polarization. Here are 

e^L) = (v,0,0,q z ) / i/Q 2 

e^x) = ( 0 , e x ) 

e n(y) = (0.4) • ( 56 ) 


The longitudinal cross section and the purely transverse unpolarized Wj = 
0.5 {e tl {x)W^' v e v {x) + e fl (y)W flI/ € l/ (y)) correspond to different combinations of W\ and 
W4, the interference term Wlt cos = —e tl (x)W^ v e l/ (L) + e ll {L)W^' v e v {x) is propor¬ 
tional to W3, and the transverse polarized term Wtt cos 2 </y- = 0.5 (e^x) W^ v e v (x) — 
e n{lj)We v (y)) is proportional to W4. 

We transform the pion variables to z, k]_, and </> 7r . With d 3 k/2E 7V = 
(n/4 k\\) dz d k\_ d</y- the five-fold differential cross section has the form 
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(57) 


(58) 


In this paper we will concentrate on <f> v integrated cross sections which depend on 
Tii and Ti 2 only. It is useful to define the combinations 


Tij = Tii and 

Hl = -Tii+ Ti 2 /2x B , (59) 

which correspond for inclusive scattering to Et = E\ and El = —E\ + E 2 l2x B . To 
select the 7-7r interaction (see Fig. 5), one must analyze the longitudinal cross section. 
Since the pion is a boson and we neglect all transverse momenta, the longitudinal cross 
section will be given entirely by the 7-7r interaction. 
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3.3. Calculation of the Structure Functions Wi(x 13 ,Q 2 ,z,k 2 1 _) 


In analogy to Eq. 53 we define a partonic tensor 

d 3 p' 


= 


1 

1 

4.X Mjy 

(2x)- 

x (g(p, 

"IP,™ 

1 

1 

4.X Mjy 

(27r) 3 

X {(({p-, 

°)\ r m 

1 

i 

2x Mjy 

(27t) 3 


£ 


(q{p XM^)l J em(0)kCp,^)) 


\4 c4/ 


£ (< l(p',<r')v(k)\J? m \q(p,<r)) 

( 7 ,( 7 ' 

(o)|?(pV'M&)) s ( p /2 - m 2 ) 

8(p' 2 -M 2 ) . 


(60) 


Here the last equation defines , which can be calculated from the Feynman dia¬ 
grams of Fig. 5 with pseudoscalar coupling. The finite size of the pion is taken into 
account by the pion electromagnetic form factor F V (Q 2 ) = m 2 /(Q 2 + m 2 ). 



+ 



u(p) d(p r ) 

(t-channel) 


Figure 5. Tree level Feynman diagrams to the process —> g?7t + . Shown are the 

contributions from the s-channel, the u-channel, and the t-channel. 


16 



In addition, we suppress all off-shell propagators of intermediate state particles with 
monopole form factors: 


g(p 2 ) 


A 2 

9 * qq A 2 + (s - Mq) 
A 2 

A 2 

9 * qq A 2 — (t — to 2 ) 


in the s-channel 
in the u-channel 
in the t-channel . 


In general, the pion-quark-quark vertex depends on three momenta, g(p 2 ,k 2 , p' 2 ). In 
the three diagrams always two of the three particles are on-shell. We simplify the 
dependence on this off-shell momentum by parametrizing it with a monopole form 
factor containing one scale A. The role of A has already been discussed in Section 2. 
The amount of a pion cloud in the u-quark depends logarithmically on A. We will use 
the value A = 3.4, which gives the best fit to the Gottfried sum rule and spin sum rule, 
see Section 2. 

The deviation of the pion form factor F v from unity could lead to a violation of 
gauge invariance. One can, however, imagine adding a contact term of the form [28] 
V 2 e g^qq Cjr u{p ', s') 75 (F V (Q 2 ) — 1) (q^/Q 2 ) u(p, s) to restore gauge invariance. Since 
e • q = 0, this term does not contribute to the calculation in Feynman gauge. A similar 
procedure has been used in Ref. [29]. 

For the process —>■ dir + we calculate M^ v = | Tr[M^(/5 + Mq) M v (ft + Mq)]. 

With q u = 2/3 and qd = — 1/3 as quark charges we get 


ikP = V2t 


9^qq 


75 + 4 + m q )y 

— Trm, — + 


+ <ld 


7^0 - jk + Mq)7, 




1 it 2 


* + F.( <V)4h2i 


( 61 ) 


The resulting tensor M^ u can be reduced according to its Lorentz structure into four 
structure functions Wq, W 2 , IF 3 , W 4 in the same way as pion production on the proton 
(cf. Eqs. 54,55). The equivalent polarization vector is T^ = x ( P^ — q^ P ■ q/q 2 ). 

The structure functions Wi are listed in the Appendix. By folding the quark struc¬ 
ture functions Wi with the quark probability distributions we obtain the hadronic 
structure functions Wi(x, Q 2 , z, k]_) 

Wr w++X = £ dx [N u (x)W? + N s (x)wf] 

P M 2 

wr v +X = l dxx2 J^ [N u (x)W-+ N d -(x)W 2 d ] 
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1 


W P^+ + X 

wr w++x 


Ax x 

o M ( 
1 


N 


N u (x)W : 3 “ + N s (x)W : 1 


0 


d.r A),(i)»7 + A'j(i)»4 


(62) 


These define the semi-inclusive structure functions (Eq. 58) and the differential cross 
section for 7 r-production (Eq. 57). 


4. The Cross Section for ep — >• e , 7r + X 


The five-fold differential cross section can be evaluated according to Eq. 57 with the 
four different structure functions Tii. For a first comparison with future experimental 
data it will be necessary - for reasons of statistics - to reduce this five-dimensional 
space to the 3 dimensions x B , Q 2 , and z by integrating out and k\\ 


1 


2-7T G_|n 


Hi {x Bl Q 2 , z) := - j dX j dk 2 L / H t (x B ,Q 2 ,zXl) 

0 0 

^ j_ |max 

= tt J dk]_ 7i t (x B} Q 2 } z } k 2 ± ) , 

0 


(63) 


where k 2 L | max is obtained from Eq. 50. The cross section is then given by 


dcr 

dx B d y dz 
dcr 

dx B d Q 2 dz 




(64) 


In Fig. 6 we show the two different contributions from Hi and H 2 to the differential 
cross section of Eq. 64 as a function of z. As parameters we choose E = 27 GeV, 
x B = 0.10, and Q 2 = 1 GeV 2 . The latter value is chosen as a compromise: Q 2 has to 
be minimal in order not to resolve the internal structure of the pion. The detection 
of the scattered electron, however, places a lower bound on Q 2 . The value of x B will 
become clear below. 

The contribution to the cross section from Hi is significantly lower than that from 
H 2 . This is due to the kinematic factor x B y 2 } which for our kinematic conditions is 
only about 0.004. Furthermore we note that the pion cloud in t-channel production 
contributes significantly to H 2} especially for large values of z. As discussed at the 
end of Section 3, this contribution is purely longitudinal owing to the pseudoscalar 
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Figure 6. Contribution from Hi and H 2 to the cross section for semi-inclusive pion 
production. The values of Q 2 and x B are fixed at 1 GeV 2 and 0.10, respectively. 

nature of the pion. This is obvious from Fig. 6, which shows a negligible contribution 
to Hi = Hx- 

Of further interest is obviously the contribution to the cross section from longitu¬ 
dinal and transverse photon polarization, respectively. Decomposing the cross section 
into longitudinal and transverse structure functions we obtain: 

= ip! {2(1 _^ L + (1 + (1 _^ lM , 

The two contributions from Hx and Hx are plotted in Fig. 7. Clearly the contribution 
from Hx at ~ —> 1 is dominated by the Cchannel, whereas the 5-channel contributes 
strongly to Hx- Thus a determination of the longitudinal part of the cross section 
and correspondingly the longitudinal structure function will allow to clearly isolate 
the presence and strength of a pion cloud. A recent calculation [25] of hadronization 
effects of the struck quark also showed evidence for a peak at large values of ~ in the 
longitudinal structure function. 

It should be stressed, that we require a large invariant mass of the hadronic system 
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excluding the pion, such that the quasi-elastic, channel 7 *p —> K + n does not contribute. 
As an example, for x B = 0.15 at E e = 27 GeV this hadronic mass is only 1 GeV 
(2.4 GeV) for ~ — > 1 (~ — > 0), with the hadronic mass increasing with decreasing x B . 
Thus a measurement at larger center-of-mass energy is desirable. 
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0.2 

: da / dx B dQ 2 dz (H T only) 
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Figure 7. Contribution from Hl and Ht to the cross section for semi-inclusive pion 
production. The values of Q 2 and x B are fixed at 1 GeV 2 and 0.10, respectively. 


Finally the question arises which value of x B is best suited to extract the signal from 
the underlying events. To this end we show in Fig. 8 the cross section as a function of 
~ for four different values of x B . Clearly the signal at large values of ~ increases with 
decreasing x B . In addition, the cross section is largest for small values of x B . Thus it 
is desirable to perform this experiment at the smallest possible value of x B . Again, for 
a typical detector the smallest value is dictated by the smallest value of Q 2 accessible 
in the experiment. Using the HERMES detector as an example, the minimum x B is 
about 0.03. 

Furthermore, the value of x B in the quark parton picture is related to x B ~ 
m(parton)/Mjv, where the structure function F- 2 has its maximum. Since in our case 
the parton is a pion we should choose x B ~ 0.15. As a compromise we use x B = 0.10. 
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z 

Figure 8. Gross section for semi-inclusive pion production as a function of the param¬ 
eter x B . The value of Q 2 is fixed at 1 GeV 2 . 


4.1. Comparison of our Model with the LEPTO Event Generator 

Experimental data on semi-inclusive pion production off nucleons is rather scarce. On# 
of the main reasons for this lack of data was the limited luminosity and limited hadron 
identification of older experiments [30]. ft is the new generation of experiments like 
the HERMES experiment [31] at HERA which have the possibility to operate at a high 
intensity circulating electron beam, which, together with a storage cell for hydrogen (or 
higher A-atoms), results in luminosities of C ~ 1 X 10 31 /cm 2 s for a polarized target. 
In case of an unpolarized target the luminosity may be increased by another factor 
of 10. Such high luminosities, together with the HERMES detector’s capability to 
measure with large efficiency hadrons produced in electron-proton collisions results in 
the possibility for a detailed study of semi-inclusive measurements like the one proposed 
in this paper. 

In the following we will compare the predictions of our model for the pion cloud of 
the quark with the results of a simulation using the Monte Carlo program LEPTO [32], 
which generates deep inelastic lepton-nucleon scattering events and hadronizes the final 
excited state using the Lund string model. We will demonstrate that our model yields 
cross sections of the same order as those simulated with LEPTO. Thus it should be 
possible to extract from data the amount of the quark’s pion cloud. Of particular help 
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will be the (near) vanishing of the longitudinal structure function, a fact known from 
several measurements. 

In Fig. 9 we compare the semi-inclusive cross sections obtained from our model 
(dashed line) and from a simulation using the fastest 7 r + from the LEPTO program 
(solid line). It is evident that the best ‘signal to noise’ at high ~ is for low values of x B . 
This is particularly true for the pion contribution from the f-channel, as can be seen 
from a comparison of the shapes of the “-distributions in Fig. 9 with those in Fig. 6. 
In addition it should be noted that the cross section is largest for small values of x B . 


da / dx B dQ 2 dz [ jub / GeV 2 ] 




Figure 9. Comparison of semi-inclusive cross sections obtained from our model (dashed 
line) and from a simulation using the fastest 7r + from the LEPTO program (solid line). 
The value of Q 2 is fixed at 1 GeV 2 . 


Also note that the LEPTO event-generator produces the expected peak at large 
~ from the elastic channel — > K + n for a value of x B = 0.15, corresponding to 
x B ~ m K /M n. In our model the wide s-channel ‘continuum’ at large x B values is 
equivalent to the hard fragmentation process, where one quark radiates a hard gluon 
which splits into a quark-antiquark pair; the antiquark then combines with the quark 
to form a pion [25]. Since for large x B values the hadronic mass excluding the pion 
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is not large compared to the nucleon mass, the calculated cross section represents an 
average cross section over the produced resonances. 

To get a feeling for the statistical accuracy of data, which can be accumulated in one 
month of data taking with an efficiency of 30% and a luminosity of C = 1 X 10 31 /cm 2 s, 
corresponding to f Cdt = 8 X 10 6 /<b _1 , we show in Fig. 10 the expected ~-spectra 
based on the LEPTO program and on our model. The data are integrated over Q 2 
from 1 to 2 GeV 2 and over x B from 0.05 to 0.15. A total of 270k events can be expected 
from LEPTO, whereas the pion model predicts 75k events. Clearly the statistical error, 
which is barely visible in Fig. 10, will allow a detailed shape comparison between data 
and LEPTO Monte Carlo, in order to ascertain a possible contribution from the pion 
cloud of the quark. 



Figure 10. The expected ~-spectra based on the LEPTO program and on our model. 
An integrated luminosity of f Cdt = 8 X 10 6 is assumed. The data are integrated 
over Q 2 from 1 to 2 GeV 2 and over x B from 0.05 to 0.15. 


However, the LEPTO Monte Carlo may not even describe the low ~ distribution 
of events, which we believe are dominated by normal fragmentation events. In this 
case any extrapolation to large ~ and any conclusion on a pion cloud contribution will 
be at least doubtful. In this case it will be necessary to separate from data directly 
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the longitudinal structure function. This requires, however, additional data taking 
at about 15 GeV, again with an integrated luminosity of f Cdt = 8 X 10 6 /rb _1 . A 
preliminary analysis [33] shows, that with two such data sets, one at 27 GeV, the other 
at 15 GeV, a determination of the longitudinal structure function will be possible with 
an accuracy of about 20% in each of 100 bins of x. This will be more than sufficient in 
order to be able to extract the prominent feature of our model for the pion cloud of the 
constituent quark: a strong peaking near x —> 1 of the longitudinal structure function. 


5. Conclusion 

Employing a model where the physical quark fluctuates with some probability to a 
quark plus a pion we have calculated the differential cross section for semi-inclusive 
pion production in electron proton reactions. As expected, the biggest contribution 
from the pion cloud to the cross section is at small values of x B ~ 0.1 and large values 
of x = E 1T /v —> 1. In particular the sole contribution to the longitudinal structure 
function will allow an experimental determination of the magnitude of the pion cloud 
of the quark. 
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6. Appendix 


In the following we list the four independent structure functions Wi for the process 
7 *u —> d'K + \ 
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Next we list the four independent structure functions Wi for the process 7 *d —> utt h 
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with the modified Mandelstam variables 
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-1/3 as quark charges. 


The factor of proportionality f3 is given by 
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